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We show that a simple, modified version of the Mag- 
netorotational Instability (MRI) can develop in the outer 
liquid core of the Earth, in the presence of a background 
shear. It requires either thermal wind, or a primary in- 
stability, such as convection, to drive a weak differential 
rotation within the core. The force balance in the Earth's 
core is very unlike classical astrophysical applications of 
the MRI (such as gaseous disks around stars). Here, the 
weak differential rotation in the Earth core yields an in- 
stability by its constructive interaction with the planet's 
much larger rotation rate. The resulting destabilising 
mechanism is just strong enough to counteract stabiliz- 
ing resistive effects, and produce growth on geophysically 
interesting timescales. We give a simple physical expla- 
nation of the instability, and show that it relies on a force 
balance appropriate to the Earth's core, known as mag- 
netostrophic balance. 



1. Introduction 

The Magnetorotational Instability (MRI) is important 
for differentially rotating astrophysical objects such as 
gaseous disks, because it forces a breakdown of laminar 
flow into turbulence, producing the enhanced dissipation 
and transfer of angular momentum necessary for accre- 
tion of matter onto the central massive object (Balbus 
and Hawley, 1991). Disk investigators began with purely 
hydrodynamical candidate mechanisms, but now the cen- 
ter of interest is squarely upon magnetohydrodynamics. 
Interest in the Earth's core, by contrast, has been mag- 
netic almost from the start. The principal problem, of 
course, has been to understand how fluid motions in 
the core generate a magnetic field. This is often ap- 
proached via kinematic processes in which a hydrody- 
namically turbulent fluid exponentially amplifies a very 
weak seed magnetic field. Here, we argue of a version of 
the MRI, dynamically coupling both the velocity and the 
magnetic field, can also provide a mechanism for linear 
instability in the Earth's core. While this mechanism is 
not meant to serve as the primary dynamo process, it 
could be an important source of secondary instabilities 
and magnetic secular variation. 



!MAG (ENS/IPGP) 
LRA, Departement de Physique 
Ecole Normale Superieure 
24, rue Lhomond 
75231 Paris Cedex 05, France. 

2 LERMA, Observatoire de Paris, CNRS/UMR8112. 

3 Institut de Physique du Globe de Paris, 
CNRS/UMR7154. 



Copyright 2009 by the American Geophysical Union. 
0094-8276/09/85.00 



1 



X - 2 



PETITDEMANGE ET AL.: MRI IN THE EARTH'S OUTER CORE 



We are hardly the first authors to study the effects 
of differential rotation on the dynamical stability of the 
gcodynamo (see e.g. Acheson, 1983; Ogden and Fearn, 
1995; Fearn et al. 1997). But in previous calculations the 
emphasis has been upon purely azimuthal fields, nonax- 
isymmetric disturbances, and magnetic instabilities. In 
this work, the dynamical focus is much different. Here 
the magnetic coupling is to the poloidal field components, 
axisymmetric disturbances are front and center, and the 
instability, while relying on the presence of a magnetic 
field, has its seat of free energy entirely in differential 
rotation. The MRI is a somewhat novel concept in this 
context, and is worthy of study in isolation. Fortunately, 
it can be understood in very direct and simple physical 
terms. 

The Earth's core is, by comparison to accretion disks, 
a relatively small object, in which resistive effects are on 
an equal footing with dynamical processes. The rotation 
properties of the Earth also significantly differ from those 
of an accretion disk. To leading order, they correspond 
to solid body rotation, with only a weak differential rota- 
tion. At first sight, it is far from obvious that the Earth's 
core is a venue for the MRI. 

The purpose of this Letter is to show that, despite its 
weak differential rotation and significant resistive effects, 
the Earth's core can in fact host the MRI. We derive 
a WKB dispersion relation relevant to this asymptotic 
regime (i.e. magnetostrophic balance), and demonstrate 
excellent agreement between this local description and 
global numerical simulations in spherical geometry. The 
final section speculates on potential applications of our 
results. 

2. Modelling 

The physical parameter regime relevant to the Earth's 
outer core dynamics is relatively well constrained. In 
order to investigate the MRI in a conducting fluid, it 
is also necessary to provide a reasonable approximation 
to the basis flow profile over which the instability may 
develop. Apart from very restricted components of the 
velocity field, imaging of the core flow is limited to sur- 
face flows (e.g. Bloxham & Jackson 1991; Amit & Olson, 
2004; Eymin & Hulot, 2005). Moreover, while secular 
magnetic field variations allow for an order-of-magnitude 
estimate for the zonal flow (Pais and Hulot 2000), it is 
not possible at present to achieve a direct, detailed or 
reconstructive imaging of the global differential rotation 
of the Earth's core. Asymptotic studies early identified 
the essential role of the zonal shear in the geodynamo 
(Taylor, 1963), and this has remained central to our un- 
derstanding of outer core dynamics. Numerical models of 
convection in rotating spheres (e.g. Christensen, 2002) or 
in convectively driven dynamos (e.g. Aubert, 2004), also 
observed the presence of a strong zonal shear. Besides 
surface flow reconstructions, the only (indirect) observa- 
tional evidence for such shear is inferred from seismolog- 
ical data that have been interpreted as rotation of the 
solid inner core at a rate of about 0.15° per year rela- 
tive to the mantle (Vidale, Dodge & Earle, 2000; see also 
Dumberry 2007). But is this change in the angular ve- 
locity profile spread across the whole core radius or, say, 
localised in a narrow shear layer? There are no observa- 
tional constraints on the way the corresponding jump in 
angular velocity is actually accommodated by the flow. 



PETITDEMANGE ET AL.: MRI IN THE EARTH'S OUTER CORE 



X- 3 



3. Magnetostrophic-MRI local description 

Consider the stability of a differentially rotating fluid, 
with angular velocity Q a function of s and z in a standard 
cylindrical coordinate system (s, (f>,z). The density p, 
pressure P, magnetic field B are functions of s and z 
in the equilibrium state, and of course depend on time t 
as well when the equilibrium is perturbed. The fluid is 
characterized by a kinematic viscosity v and resistivity r\. 
We work in the WKB limit, assuming that all linearized 
perturbations have the space-time dependence 

Q oc exp(ifc s s + ik z z + at) (1) 

where Q is an infinitesimal (Eulerian) disturbance in a 
fluid quantity. Our starting point is the general disper- 
sion relation of a constant density fluid adapted from 
Menou, Balbus, & Spruit (2004): 



£.2 <J V' / °V : 



where 



and 
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4fi 2 (fc- Va) 2 =0, (2) 



Sv = [{a + f)k 2 ){cr + uk 2 ) + (k ■ Va) 2 } , (4) 

a 2 v = [(a + V k 2 ) 2 + (k- Va) 2 ] . (5) 
The Alfven velocity Va is defined as 

Va = -%= (6) 



tip 

where \i is the magnetic permitivity. Finally, we note 
that is the full angular velocity, i.e. if flo denotes the 
Earth rotation rate, and V the velocity in the rotating 
frame, then Q = Qo + V$/s . 

To see whether the unstable MRI modes can be present 
in a very simple model, we consider the case of an ini- 
tial poloidal magnetic field which is locally vertical and 
Q is a function only of s. (We stress here that we use 
a vertical magnetic field for the sake of simplicity of the 
presentation, but that the overall conclusion can be eas- 
ily extended to the case of more general applied fields.) 
Then 

2 d£l 2 \ 2 



v^n 2 ) = - \4£l 2 + s^J = -k. (7) 

For the current application, the shear gradient is very 
small compared with 4Q 2 , and k, which is known as the 
"epicyclic frequency" in the astrophysical literature, is 
very nearly 2fi. With these assumptions, the dispersion 
relation becomes 

% Ua + V k 2 )(a + uk 2 ) + (k ■ Va) 2 } 2 + 4Q 2 (a + r,k 2 ) 2 + (k ■ V A fs^- = . 
ki 1 J as 

(8) 

This is the form that we have used (see below) for com- 
parison with numerical simulations. 

In geodynamo applications, the viscosity is very small 
compared with the resistivity (the ratio is 10 -8 — 10 -6 ) 
and it may be dropped. The instability of interest arises 



X - 4 



PETITDEMANGE ET AL.: MRI IN THE EARTH'S OUTER CORE 



from the final term of this dispersion relation, if the an- 
gular velocity decreases outwards. This is a small term 
in the sense that fi is regarded as large, and we expect 
that the growth rate a will itself be very small compared 
with fi. We may therefore drop all a terms from our dis- 
persion relation, except for those which are multiplied by 
Q, 2 . (This procedure may be formalized by normalizing 
a with respect to \sdQ,/ds\ and then treating 

Q, ds 

as a vanishingly small expansion parameter.) We obtain 
^(k ■ Va) 4 + 4{l 2 (a + n k 2 ) 2 + (k ■ V A fs ^f- = 0(9) 



This relatively simple dispersion relation has a cor- 
respondingly simple physical interpretation. Consider 
magnetostrophic balance with an axial wave number 
(k = fee*), 

2fi x v = -vn + — (B-V)B, (10) 
MP 

where n includes the material and magnetic pressures 
as well as the centrifugal potential, along with the full 
induction equation 

(|^ + F-v)s = (B-V)y + J? V 2 B. (11) 

The equilibrium profile is to leading order V4, = s[Qo + 
(s — so) dQ/ds], and B — Be z . Consider linear velocity 
perturbations v s and v^>, and magnetic field perturba- 
tions b s and b# of the form exp(crt + ikz). In magne- 
tostrophic balance, the radial and azimuthal linearized 
equations are respectively 

-2fiv* = ^ 6 S , (12) 
pop 

2Vv s = ^ & . (13) 
Pop 

The same components of the induction equation are 

(a + rik 2 )b s =ikB v s , (14) 

(<r + r/k 2 ) b<p = ikBo v<p + s^-bs . (15) 
Combining (14) and (13) leads to 

(a + r;fc 2 )6 s = -^(fcVA)V (16) 
whereas (15) and (12) imply 

'(kV A ) 2 dO. 



(a + rjk )b<j, 



20, +S 'ds 



(17) 



Notice the absence of pressure-like perturbations for axial 
wavenumbers. The last two equations combine to give 

dQ 2 

(Wa) 4 + 40 2 (a + V k 2 ) 2 + (kV A ) 2 s^- = 0, (18) 



which is just equation (9) for axial wavenumbers. 

By restricting the perturbations exclusively to the 
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magnetic field components, equations (16) and (17) pro- 
duce a clear picture of the instability with no "out-of- 
phase" terms. Consider an outward radial displacement 
and its associated radial magnetic field. If UVa is not 
too large, the radial field is significantly sheared by the 
differential rotation to produce a negative azimuthal elon- 
gation of the field line following equation (17). This 
results in a positive azimuthal magnetic tension force, 
which must be balanced in the magnetostrophic regime 
with an opposing Coriolis force, hence with a yet greater 
radially outward displacement and magnetic field (16). 
An instability is at hand. This mechanism is illustrated 
in Figure 1, and seen through a direct simulation in Fig- 
ure 2. 

The growing solution for a is 



|fc • Va\ 

20, 



dO z 
ds 



k 



1 1/2 



k2 (k-v A y 



rjk 2 .(19) 



We wish to find the maximum growth rate of the insta- 
bility, which will be associated with a particularly wave 
vector k. The easiest way to proceed is to use the vari- 
ables 



X = k- Va, 



kV 



(20) 



with the understanding that X > 0, and the parameters 



Vj _ gg 
2r/O 2fj, prjO ' 



ds 



(21) 



A is the Elsasser number, which is typically of order unity 
for the Earth's core. Then 



20a = X(a - X 2 Y) 1 ' 2 - X 2 Y/A . 



(22) 



This is a decreasing function of Y everywhere, with a 
maximum at the lower boundary Y = 1. This means 
that the most rapidly growing wavenumber k z lies along 
the rotation axis. Thus 



2n<x (max) = X(a - X 2 ) 1/2 - X 2 /A 



(23) 



Requiring the partial derivative d/dX to vanish leads to 
the polynomial 



X 4 - aX 2 + 



= 0. 



4(1 + A" 2 ) 

The physically meaningful wavenumber solution is 
dO 



(24) 



X 



(k z v A y 



ds 



[l-(l + A 2 )- 1/2 ] , (25) 



and corresponds to a growth rate of 



dtt 
ds 



A/2 



l + yiTA 2 " 



(26) 



It is of interest to examine the range of k z in which 
instability exists. For the case k — k z , the dispersion 
relation (9) admits growing solutions for values of k 2 z less 
than 



\sdn 2 /ds\ 



V 2 +4Q 2 ri 2 /V 2 



(27) 
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Notice that when Va is very large, k z is very small be- 
cause large k z perturbations are stabilized by magnetic 
tension forces. Conversely, when Va is tiny, resistiv- 
ity damps large wavenumber perturbations. The max- 
imum value of k z allowing the greatest range of unsta- 
ble wavenumbers corresponds to A = 1. Hence, Elsasser 
numbers of order unity naturally emerge in an MRI- 
influenced dynamo. 

At this point it is helpful to have some explicit num- 
bers. We take fi = 7.27 x 1(T 5 s _1 and sdfl/ds ~ 1(T 10 
s _1 , the latter corresponding to a lower bound for the an- 
gular velocity gradient in which the characteristic length 
scale associated with shear in the core is the outer core 
radius itself (worse case MRI scenario) . With A = 1 and 
r\ ~ 1 m 2 s~ , the Alfven velocity is 0.012 m s _1 . From 
equation (25), we find 

W A ^ 4.6 x 10~ 8 s _1 

which translates to a wavelength of some 1600 km. The 
characteristic growth time from equation (26) is then 
about 1500 years. If, on the other hand, we assume that 
a similar jump in angular velocity is accomodated accross 
a narrow shear layer (such as the Stewartson layer, Stew- 
artson 1966), then the characteristic growthtime could be 
as short as a year. 

What is the relationship between the unstable mode 
considered here and the maintenance of the classical Tay- 
lor constraint? It is well known that in the magnc- 
tostrophic limit, the geostrophic component of the flow 
must adjust to ensure Taylor's constraint (Taylor, 1964). 
A modification of the non-geostrophic flow similar to the 
one envisioned here could in principle rapidly alter the 
zonal shear on which the instability itself relies. How- 
ever, in the particular configuration investigated, the 
wavenumber is axial. Thus, at least in the linear phase 
of the disturbance, the Taylor constraint remains unaf- 
fected. 



4. Direct numerical simulations 

To demonstrate numerically the above mechanism, we 
base our study on a very idealized model of the Earth 
core, chosen for illustrative purposes. We consider sim- 
ple spherical Couette flow driven by enforcing differential 
rotation between the inner core and the mantle. The ap- 
plied magnetic field Bo is vertical and uniform. In the 
simpler hydrodynamical case, a strong shear layer will 
develop on the cylinder tangent to the inner core (Proud- 
man, 1956; Stewartson, 1966). With the particular choice 
of a vertical applied field, this shear is only slightly mod- 
ified by MHD effects (see however Dormy et al., 1998). 
In what follows, we will work in time units of (2f2o) _ 
and space units of r , the radius of the outer sphere. It is 
also convenient to introduce the dimensionless velocity 

Uo = ^r- , (28) 

where Vo is the unperturbed velocity measured in the 
frame rotating at £7o- As before, the linear perturbed 
velocity is v. For the perturbed velocity v and magnetic 
field b, we introduce dimensionless variables v' and b': 



v' = v/(2r o Q ) , b' = b/B . (29) 
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The governing fluid equations may then be written 

^ + (U ■ V)v' + (v' ■ V)U = -Vtt + EAv' + |^-(V x b') x e z - e z x v' , 

(30) 

^=Vx(v' xe z + U xb') + -^-Ab' , (31) 
at Pm 

V • v' = V • b' = . (32) 

In addition to the Elsasser number A, we have intro- 
duced the standard Ekman and magnetic Prandtl num- 
bers, which are respectively: 

E = u/(2n r 2 o ) Pm = v/r), (33) 

and 7T is the perturbed value of the dimesionless pressure. 

As a model for the Earth's core, we investigate a spher- 
ical shell. The angular velocity Q.q is defined by the outer 
sphere. In this reference frame, the inner sphere rotates 
with angular velocity Qi > 0, in order for the differen- 
tial rotation to decrease outward. We define the Rossby 
number as 

Ro = Q z /Qq. (34) 

Let us now compare the local WKB analysis (8) with 
the global growth rate obtained numerically. We first 
need to compute the dimensionless steady velocity pro- 
file Uo- If the magnetic field is weak enough, this can 
be obtained as a solution of the hydrodynamic problem. 
The results are displayed in Table 1. If the Elsasser num- 
ber (21) becomes large, the field obviously affects the 
steady solution, and we then perform simulations with 
the steady MHD state as initial configuration, see Ta- 
ble 2. 

To best illustrate the above theoretical analysis, we 
use numerical values as close as feasible to the regime 
described. Simulations use Pm — 0.5, the differen- 
tial rotation is weak compared to the planet rotation 
(Ro « 1), and the Ekman number is decreased to small 
values. To accommodate small Ekman numbers, we have 
implemented an axisymetric version of the Parody code 
(Dormy et al., 1998, Christensen el al., 2001, and later 
collaborative developments) with radial resolution rang- 
ing from 500 to 1000 points in radius and between 200 
and 500 harmonics. 

Our results are summarized in Table 1. In agree- 
ment with theoretical expectations, the radial wavevector 
k. 9 varies approximately in proportion to the dominant 
shear, i.e. E 1 ^ 4 . The value of k s is taken from numeri- 
cal simulations by computing the width of the unstable 
mode at mid-intensity, and the value of k z is obtained 
by direct fourier transform. It is found that the most 
unstable mode's half wavelength generally occupies the 
full radial extent of the sheared region. The inner core 
and the associated equatorial singularity of the Ekman 
layer are the natural modal boundaries; the equatorial 
singularity divides the tangent cylinder into independent 
domains. Finally, to limit the stabilizing effect of diffu- 
sion, we must use E much less than unity. 

The induced azimuthal magnetic field clearly repro- 
duces the mechanism of the magnetostrophic MRI (see 
figure 1), despite the presence of finite inertial and diffu- 
sive effects (which were neglected in the theoretical calcu- 
lation) , and the complications associated with a bounded 
spherical domain (also neglected). 

Some of these cases involve large values of the Elsasser 
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number, one may then worry that they are not fully self 
consistent. To check this, we show, in Table 2, the in- 
stability parameters of such strong field profiles, which 
now also have a (weak) dependence on A. For the field 
geometry used here, however, both the rotation profile 
and the applied field require relatively little adjustment. 
Agreement with the local description is also obtained for 
such configurations. 

5. Conclusions and Discussion 

At the very least, the magnetostrophic MRI discussed 
here seems to be an efficient regulatory mechanism in 
the outer core, providing a back reaction against the 
build up of angular velocity at smaller radii. This could 
be tested in full scale numerical simulations of spherical 
shells whose gross aspect and parameters could be regu- 
lated either to admit the MRI or not. 

More interestingly, might there be traces of the mag- 
netostrophic MRI in geomagnetic observations? Because 
of the severe limitations in our knowledge of differential 
rotation within the Earth's core, it is not possible to offer 
any kind of definitive treatment. It is, however, interest- 
ing to note that some axial variation in the zonal flow is 
required to account for changes in the length of the a day 
over millenial timescales (Dumberry & Bloxham 2006). 
In that respect It is interesting that axial wavenumber 
disturbances grow particularly rapidly in the model stud- 
ied here. 

Another temptative application is based on the obser- 
vation that the growth rate in question is directly propor- 
tional to the shear, which could be highly localized. Thus, 
growth times considerably shorter than 1500 years are 
possible: in principle, shear times of order years could be 
produced. This suggests a possible connection with the 
observations of rapid events of internal origin, known as 
"geomagnetic jerks" or "geomagnetic impulses". These 
rapid changes in the secular variation were detected and 
characterized from observatory data (Courtillot et al, 
1978; Le Mouel et al, 1982) and found to correspond to 
localised patches of rapid field variation at the surface of 
the core, using global field models (Dormy and Mandea, 
2005). We reemphasize, however, that the establishment 
of a compelling link between MRI unstable configurations 
in the Earth core and geomagnetic impulses will require 
considerably more effort. Nevertheless, the MRI offers for 
consideration a novel magnetohydrodynamic mechanism 
that is able to induce significant field variations over a 
range of observationally interesting timescales. 
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Table 1. Results obtained for the most unstable mode in 

a direct integration of the governing equations. We list the 

growth rate (cTnum) and the wavenumber (k s and k z ). We 

also list the growth rate (<Tth ) obtained directly from the local 

dispersion relation (8) with these wavenumber values. The 

variables E, Ro, and A are defined in equations (21), (33) and 

(34). The numerical agreement is excellent. 
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10 


7.1 x 10" 


-3 


7.6 x 10~ 


-3 


vr/0.11 


tt/0.1 






5 


6.0 x 10" 


-3 


6.4 x IO" 


-3 


tt/0.1 


tt/0.071 






2 


3.0 x 10" 


-3 


4.8 x 10- 


-3 


tt/0.1 


tt/0.067 






1 


2.24 x 10 


-3 


3.0 x 10- 


-3 


tt/0.1 


tt/0.058 






0.5 


1.16 x 10 


-3 


1.36 x 10 


-3 


tt/0.11 


tt/0.11 




0.005 


0.3 


6.96 x 10 


-4 


2.03 x 10 


-4 


tt/0.1 


tt/0.14 




0.0025 


25 


2.0 x 10" 


-3 


1.93 x 10 


-3 


tt/0.1 


tt/0.25 




0.0075 


25 


6.2 x 10~ 


-3 


6.6 x 10- 


-3 


tt/0.1 


tt/0.2 


2.5 x icr 7 


0.0025 


25 


2.0 x 10" 


-3 


3.2 x 10- 


-3 


tt/0.08 


tt/0.15 


io- 7 


0.0025 


1 


3.0 x 10~ 


-4 


2.2 x 10- 


-4 


tt/0.06 


tt/0.04 






25 


4.0 x 10" 


-3 


4.8 x 10- 


-3 


tt/0.06 


tt/0.083 



Table 2. Results obtained using the steady MHD state as 
initial configuration. Computations are performed with Ro = 
5 X 10~ 3 , Prn = 0.5 and E = 2.5x 10~ 7 . As for Table 1, both 
the numerical parameters of the most unstable mode and the 
corresponding theoretical growth rate are reported. 



2A 


fth 






k z 


3 


5.06 x 10" 3 


5.70 x 10" 3 


tt/0.92 


tt/0.056 


1 


2.06 x 10~ 3 


1.84 x 10~ 3 


tt/0.95 


tt/0.056 


0.3 


8.44 x 10~ 4 


2.76 x 10~ 4 


tt/0.10 


tt/0.11 
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Figure 1. Development of the magnetostrophic MRI. (a) Unperturbed field line, (b) Field line is distorted by a radially 
outward displacement, and subject to velocity shear, (c) Field line develops azimuthal tension which is immediately 
compensated by the Coriolis force. This compensating force requires a further displacement in the same sense of the initial 
outward radial extension, and the instability proceeds. 





Figure 2. Snapshot of the most unstable mode in a direct simulation with E = 5 x 10~ 7 , A = 3 and Pm 
The physical mechanism of the magnetostrophic-MRI can be traced in the phase shifts between v 3 , b 
quarter-period ahead of v 3 , while is exactly out of phase by half a period with b s . 
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and b^: b s is a 



